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Resumen 


Cortés-Martínez, F., Treviño-Cansino, A., Espinoza-Fraire, A. 
T.,, Sáenz-López, A., Alcorta-García, M. A., González-Barrios, 
J. L,, Martínez-Rojas, R., € De la-Cruz-Acosta, EF. J. (julio- 
agosto, 2017). Optimización en el diseño de un sistema de 
tratamiento de aguas residuales integrado por tres lagunas 
de estabilización. Tecnología y Ciencias del Agua, 8(4), 139-155. 


El artículo presenta un modelo matemático para la 
optimización de un sistema lagunar integrado por tres 
lagunas: anaerobia, facultativa y de maduración. El 
modelo matemático de optimización tiene como objetivo 
la minimización del área ocupada por el sistema lagunar, 
y presenta como restricciones la calidad del agua en el 
efluente: demanda bioquímica de oxígeno y el número de 
coliformes fecales; además, se restringe la dispersión (mayor 
o igual a 0.05), con el propósito de favorecer el flujo disperso 
en la laguna facultativa y en la laguna de maduración. Las 
variables cambiantes en el modelo de optimización son las 
siguientes: tiempo de retención y número de mamparas 
tanto de la laguna facultativa como la de maduración, y la 
relación largo ancho para la laguna anaerobia. Se diseñó un 
sistema de lagunas, luego se aplicó el modelo matemático 
propuesto. Los resultados muestran una disminución 
del área de 15.16%. Los dos estudios cumplen con límites 
máximos permisibles de descarga de acuerdo con la NOM- 
SEMARNAT-1996: calidad del agua residual tratada para el 
vertido a los cuerpos receptores. 


Palabras clave: dispersión, optimización, restricciones, 
sistema lagunar y tiempo de retención. 


Abstract 


Cortés-Martínez, F., Treviño-Cansino, A., Espinoza-Fraire, A. 
T., Sáenz-López, A., Alcorta-García, M. A., González-Barrios, ]. 
L., Martínez-Rojas, R., €' De la-Cruz-Acosta, E J. (July-August, 
2017). Optimization in the design of a wastewater treatment system 
integrated by three stabilization lagoon. Water Technology and 
Sciences (in Spanish), 8(4), 139-155. 


Article presents a mathematical model for the optimization of a 
lagoon system composed of three lagoons: anaerobic, facultative 
and maturation. The mathematical model of optimization aims at 
minimizing the area occupied by the lagoon system, and has as 
restrictions the water quality in the effluent: biochemical oxygen 
demand and the number of fecal coliforms; In addition, the dispersion 
(greater than or equal to 0.05) is restricted, in order to favor the 
dispersed flow in the facultative lagoon and the ripening lagoon. 
The changing variables in the optimization model are the following: 
retention time and number ofscreens, both facultative and maturation 
lagoons, and the long wide ratio for the anaerobic lagoon. A lagoon 
system design was performed, then the proposed mathematical model 
was applied. The results show a decrease of the area of 15.16 percent. 
The two studies comply with maximum permissible discharge limits, 
according to NOM-SEMARNAT-1996: quality of wastewater 
treated for discharge to receiving bodies. 


Keywords: dispersion, optimization, constraints, lagoon system 
and retention time. 
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Introducción 


Una laguna de estabilización es un reservorio 
de agua con poca profundidad. Este tipo de 
tratamiento representa una ventaja en relación 
con los demás sistemas de tratamiento: cons- 
trucción sencilla, y fácil operación y manteni- 
miento (Metcalf €: Eddy Inc., 1991). De acuerdo 
con la Comisión Nacional del Agua (Conagua), 
las lagunas de estabilización se clasifican en 
anaerobias, facultativas y de pulimiento, o de 
maduración; su objetivo principal es la elimina- 
ción de contaminantes: demanda bioquímica de 
oxígeno (DBO) y coliformes fecales (Conagua- 
IMTA, 20076). Las lagunas anaerobias funcionan 
en ausencia de oxígeno y las profundidades 
varían de 3 a 5 metros. El tratamiento se lleva a 
cabo por medio de bacterias que forman ácidos 
orgánicos y bacterias metanogénicas. 

Las lagunas facultativas incluyen tres partes: 
la primera se conoce como zona superficial, en 
donde algas y bacterias tienen una relación sim- 
biótica; la zona de fondo es anaerobia, en donde 
se acumulan los lodos y son descompuestos; en 
la zona intermedia, conocida como facultativa, 
se descompone la DBO o materia orgánica, a 
través de las bacterias aerobias, anaerobias y 
facultativas; según Rolim (2000), las profundi- 
dades varían de 1.5 a 2.5 metros. 

Las lagunas de maduración o de pulimiento 
son aerobias y se utilizan al final del proceso; 
su propósito es incrementar la calidad del agua 
residual reduciendo los organismos patógenos; 
la eliminación de las bacterias depende de 
las condiciones climatológicas y ambientales, 
por ejemplo: un potencial de hidrógeno (pH) 
alto y la exposición al sol. Las profundidades 
sugeridas son de 0.6 a 1.5 metros (Rolim, 2000; 
Conagua-IMTA, 2007a). 

Los patrones de flujo dentro de las lagunas 
de estabilización juegan un papel importante: 
el sistema biológico se basa en el sistema 
hidráulico, por lo tanto es prudente llevar a 
cabo una revisión cuidadosa de la hidráulica 
de las lagunas. De acuerdo con Rolim (2000), 
los patrones de flujo y mezcla en lagunas de 
estabilización se clasifican en: a) flujo pistón, 
el caudal se desaloja en el mismo orden que 
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entró y no se establece mezcla ni dispersión; 
b) mezcla completa: los elementos del caudal 
son mezclados y su contendido es homogéneo 
en toda la laguna; c) flujo disperso, el caudal 
tiene un tiempo de retención diferente para cada 
periodo, este patrón también es conocido como 
flujo arbitrario y se localiza entre dos límites: 
flujo pistón y mezcla completa. 

Wehmer y Wilhelm (1956) llevaron a cabo 
un estudio donde relacionan la remoción de 
los contaminantes con la dispersión dentro de 
las lagunas de oxidación: ubicaron un submo- 
delo hidráulico de un reactor entre los límites 
de mezcla completa y flujo pistón. Más tarde, 
Thirumurthi (1969) analizó que en las lagunas 
no se implementaba el flujo pistón ni mezcla 
completa y que, efectivamente, sí se establecía 
el flujo disperso. Tiempo después, Rolim (2000), 
Mara (2004), y Shilton y Bailey (2006) reportaron 
que la geometría de las lagunas tenía influencia 
importante sobre el desempeño hidráulico y la 
eliminación de contaminantes; por ejemplo, una 
relación largo ancho de tres favorece el estable- 
cimiento del flujo pistón en el estanque; es decir, 
se refiere a lagunas rectangulares. Teóricamente, 
el patrón de mezcla completa se favorece cuan- 
do las lagunas son cuadradas o circulares. El 
flujo disperso está comprendido entre los dos 
patrones anteriores. 

En la figura 1 se muestra una gráfica donde 
se relaciona el coeficiente de reducción bacteria- 
na y la materia orgánica sobrante; se incluyen 
diferentes escenarios para flujo disperso (Weh- 
ner £ Wilhelm, 1956; Rojas $: León, 1990). 


Mamparas o deflectores 


Shilton y Mara (2005), así como Abbas, Nasr y 
Seif (2006) recomendaron que en el diseño se 
incluyeran mamparas o canales, pues mejoran la 
hidráulica dentro de los estanques. Los siguien- 
tes investigadores: Killani y Ogunrombi (1984); 
Pedahzur, Nasser, Dor, Fattal y Shuval (1993); 
Muttamara y Puetpaiboon (1997); Zanotelli, 
Medri, Belli-Filho, Perdomo y Costa (2002); Shilton 
y Harrison (2003); Sperling, Chernicharo, Soares 
y Zerbini (2003); Shilton y Mara (2005); Abbas et 
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Porcentaje remanente de C/Co 


Flujo pistón: d = 0 
cs ¿[| = 0,05 


sm] = 0,1 
cs [| = 0,2 
ms [ = 0,3 
d=0.4 
ms  = 0.5 
d=0.75 
mm  =1.0 
md =1.5 
—d=2.0 
n= d=3.0 
———d=5.0 


——¿-80 
===" Flujo completamente mezclado: d = 00 


Valores de K.t 


Figura 1. Reducción de coliformes fecales con flujo disperso (Wehner éz Wilhelm, 1956; Rojas €: León, 1990). 


al. (2006), llevaron a cabo estudios experimen- 
tales con deflectores. Concluyeron que con una 
longitud de 70% de los deflectores en la laguna 
se obtienen mejores resultados en el tratamiento 
del agua residual; es decir, se incrementa la efi- 
ciencia en la remoción de contaminantes. 


Acerca de los modelos matemáticos 


De la revisión bibliográfica acerca del tema se 
tiene que Killani y Ogunrombi (1984) sugirieron 
la aplicación de la investigación de operaciones 
en sistemas lagunares. Como resultado de la op- 
timización obtuvieron el costo óptimo. Luego, 
Fonseca y Fleming (1993, 1995); Oke y Otun 
(2001); Bracho, Lloyd y Aldana (2006); Winfrey 
Stronsnider, Nairn y Strevett (2010); Olukanni 
y Ducoste (2011), emplearon la investigación 
de operaciones sugerida por Killani y Ogun- 
rombi; determinaron que era posible reducir el 
costo del sistema, al igual que incrementar la 
eficiencia en el tratamiento del agua residual. 
Enseguida, Sah, Rousseau y Hooijmans (2012) 
llevaron a cabo comparaciones de modelos ma- 
temáticos para optimización y recomendaron la 
necesidad de adecuar los modelos de tal forma 


que se incluyera un mayor número de variables. 
Después, Lee y Cheong (2014) reportaron un 
análisis numérico de la hidráulica en una laguna 
de estabilización, para optimizar las condiciones 
del flujo. Los resultados: una profundidad de 
dos metros y una relación largo-ancho igual 
a cinco. Luego, Cortés-Martínez, “Treviño- 
Cansino, Luévanos-Rojas y Luévanos-Rojas 
(2013); Cortés-Martínez, Treviño-Cansino, Lué- 
vanos-Rojas, Luévanos-Rojas y Uranga (2014a); 
Cortés-Martínez, Treviño-Cansino, Kalashnikov 
y Luévanos-Rojas (2014b); Cortés-Martínez, 
Treviño-Cansino, Alcorta-García, Sáenz-López 
y González-Barrios (2015), y Cortés-Martínez, 
Treviño-Cansino, Sáenz-López, González- 
Barrios y De la-Cruz-Acosta (2016) publicaron 
modelos matemáticos para sistemas lagunares 
con diferentes configuraciones, que consideran 
como función objetivo el costo del sistema de 
tratamiento, restringido siempre por la norma 
de calidad del agua residual. 


Algoritmo para la optimización 


Se aplicará el sistema Solver de Excel, que inclu- 
ye el algoritmo gradiente reducido generalizado 
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(GRG): resuelve problemas de programación 
lineales o no lineales, y admite expresiones 
matemáticas de igualdad o desigualdad. En 
términos generales, el Solver comienza la opti- 
mización con un valor conocido (una solución), 
luego el algoritmo la mejora hasta donde le 
permiten las restricciones previamente definidas 
(Muramatsu, 2011). 

El presente documento tiene como propósito 
minimizar el área en el diseño de un sistema 
lagunar, considerando: a) una función objetivo y 
restricciones para la optimización de un sistema 
compuesto por tres reactores: anaerobio, facul- 
tativo y maduración; b) diseñar el sistema de 
estanques con el criterio tradicional; c) emplear 
el modelo matemático utilizando el sistema 
Solver de Excel, y d) comparar los resultados 
entre ambos análisis. 

El documento presenta como aportación 
un nuevo modelo matemático, que incluye la 
dispersión en las restricciones, con el propósito 
de mantener el flujo disperso. Lo anterior, con- 
siderando un sistema de tratamiento integrado 
por tres lagunas. 


Metodología 


En la primera etapa se consideró la nomen- 
clatura y descripción del método adoptado 
por la Comisión Nacional del Agua. En este 
apartado se construyó el modelo matemático de 
minimización del área ocupada por un sistema: 
laguna anaerobia, facultativa y de maduración. 
El modelo se compone de una función objetivo, 
la cual se forma por la suma de las áreas ocupa- 
das por cada una de las lagunas, y tiene como 
restricciones el número de coliformes fecales, la 
DBO en el efluente, así como la dispersión y el 
tiempo de retención. Posteriormente se analiza 
un caso, calculando el sistema lagunar con el 
método tradicional de diseño de la Conagua y 
otro aplicando el modelo matemático, donde se 
toman como variables el tiempo de retención 
y el número de mamparas tanto de la laguna 
facultativa como la de maduración, y la rela- 
ción largo ancho para la laguna anaerobia. Lo 
anterior con el propósito de encontrar el arreglo 
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más eficiente en uso de área que cumpla con las 
restricciones del modelo. Por último se discuten 
los resultados y se hace un análisis de sensibili- 
dad del modelo matemático. 

El método para el diseño del sistema de 
lagunas formado por tres reactores fue consul- 
tado en el manual Paquetes tecnológicos para el 
tratamiento de excretas y aguas residuales en comu- 
nidades rurales (Conagua-IMTA, 2007a). 


Nomenclatura 


> 
SS 
Il 


área de la laguna anaerobia (m?). 


a, = constante adimensional laguna facul- 
tativa. 

A, = área de la laguna facultativa en (m?). 

dle = área de la laguna de maduración (m2). 

a, = constante adimensional, laguna de 
maduración. 

Asu)  = área superficial de la laguna (m?). 

A, = área total (m?). 

Bs,p  = ancho superior de la laguna (m). 


prom = ancho promedio de la laguna en (m). 

CO  = carga orgánica (kg/d). 

DBOea = concentración de la DBO, en el efluen- 
te de la laguna anaerobia (mg/l). 

DBOf£f = concentración de la DBO, en el efluen- 
te de la laguna facultativa (mg/l). 

DBOem =concentración de la DBO, en el efluen- 
te de la laguna de maduración (mg/1). 

DBOia = concentración de la demanda bioquí- 
mica de oxígeno en el influente de la 
laguna anaerobia (mg/1). 

DBOf[f = concentración de la demanda bioquí- 
mica de oxígeno en el influente de la 
laguna facultativa (mg/l). 

d, = factor de dispersión de la laguna 
facultativa adimensional. 

d = factor de dispersión de la laguna de 

maduración adimensional. 

2.7182818. 


a 
Il 


ev = evaporación (mm /d). 

K, = coeficiente de reducción bacteriana 
(d*). 

Kf = constante de decaimiento a una 


temperatura en cualquiera día con 
unidades (d”). 
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L. = concentración de la demanda bioquí- 
mica de oxígeno en el influente de la 
laguna anaerobia (mg/l). 

Kta  =constante de decaimiento para la 
laguna anaerobia. 

L, = concentración de DBO, en el influente 
de la laguna facultativa (mg/l). 

L = longitud promedio de la laguna (m). 

= largo superior de la laguna (m). 

Nea  =coliformes fecales corregidos por 
evaporación en laguna anaerobia 
(NMP /100 ml). 

Nef  =coliformes fecales corregidos por 
evaporación en laguna facultativa 
(NMP /100 ml). 

Nem  = coliformes fecales corregidos por 
evaporación en laguna de maduración 
(NMP /100 ml). 

Nf/No = número de coliformes fecales en el 
efluente (NMP /100 ml). 

Ni = coliformes fecales en el influente de la 
laguna (NMP/100 ml). 

mampe Número de mamparas en la laguna 
facultativa. 

O = tiempo medio de retención hidráulico 
de la laguna anaerobia en días. 

O, = caudal en el efluente de la laguna 
anaerobia (m? / d). 

O...  = caudal en el efluente de la laguna de 
maduración (m?/d). 

O, = caudal en el influente de la laguna 
anaerobia (m? / d). 

O, = caudal en el influente de la laguna 
facultativa (m?/d). 

O,,, = caudal en el influente de la laguna de 
maduración (m?/d). 

Omed = caudal en el influente en (m? / día). 


Of = tiempo medio de retención hidráulico, 
laguna facultativa (d). 

Om =tiempo de retención en laguna de 
maduración en días. 

T = temperatura mínima media mensual 
del aire (*C). 

Talud = relación de inclinación de bordos. 

Va = volumen de la laguna anaerobia (m?). 

Vf = volumen de la laguna facultativa (m?). 


Vm = volumen de la laguna de maduración 
(m3). 
Xf = relación entre longitud y anchura de 


la laguna facultativa. 
1000  = factor de conversión. 


Av = carga orgánica volumétrica (gDBO,/ 
m*-d). 

As = carga superficial de diseño (kg /ha-d). 

X, = relación entre longitud y anchura de 
la laguna anaerobia. 

Za = profundidad de la laguna anaerobia 
(m). 

2; = profundidad de la laguna facultativa 
(m). 

Z. = profundidad de la laguna de madura- 
ción (m). 


Laguna anaerobia 


La carga orgánica se calcula mediante la 
ecuación (1). 


Carga orgánica: 
CO = Qia(DBOia)/1 000 (1) 


Carga superficial de diseño: 


Mo =20(T) -100 (2) 
Remoción de la DBO: 
% DBO removido=2T+20 (3) 


Volumen de la laguna: 
Va =LiaQia / Mo (4) 


Área de la laguna. Profundidad de la laguna 
(Z) de2a4m: 


Aa=Va/Za (5) 
Tiempo medio de retención hidráulico: 


Oa = Va/Qia (6) 
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Concentración de la DBO en el efluente de la 
laguna: 


DBO€ =(100 -%DBO removida(DBOia) (7) 
Gasto en el efluente: 
Qea = Qia —-0.001 Aaev (8) 


Remoción de coliformes fecales. Constante 
global de decaimiento: 


Kta.(d)=2.6(1.19)72 (9) 


Coliformes fecales en el efluente de la laguna: 


Ni 


Nea = (10) 
1+Kta0a 
DBO corregida por evaporación: 
DBO ea =(DBO¡a)Qiía) / O ea (11) 


Coliformes fecales corregidos por evaporación: 


Ne =(Ne)NQia) / Qea (12) 
Dimensionamiento de la laguna. Relación largo 
ancho X =2: 

Bprom = pa (13) 
Xa 
Lprom = (14) 
Bprom 
Bsup = Bprom +(Za)(Talud) (15) 
Lsup = Lprom+(Za)lXTalud) (16) 
Cálculo del área superficial: 
Asup =(Bsup)lLsup) (17) 


Diseño de laguna facultativa (flujo 
disperso: método de Yánez) 


a) Carga orgánica. Se aplica la ecuación (1) de 
la laguna anaerobia. 
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b) Carga superficial de diseño: 
As =250(1.085)" (18) 
c) Área de la laguna facultativa. 


Se calcula con base en la carga superficial me- 
diante la siguiente ecuación: 


Af =10L,Qmed / hs (19) 
d) Volumen de la laguna: 


Vf =(APZF) (20) 


e) Tiempo de retención hidráulico. Aplicar la 
ecuación (6) de la laguna anaerobia. 

f) Dimensionamiento. Relación largo ancho 
X =3. Para determinar el ancho y longitud 
promedio del estanque facultativo emplear 
las ecuaciones (13) y (14); para el ancho y 
largo superiores (B,,, y L;,,), considerar las 
fórmulas (15) y (16); la inclinación de los 
taludes 2:1 y la profundidad Z =1.5 m. 

g) Cálculo del área superficial. Aplicar la 
expresión (17). 

h) Gasto en el efluente. Aplicar la expresión (8) 

i) Remoción de coliformes fecales. La ecuación 
considera la relación largo ancho (X). 
Incluyendo deflectores con una longitud 
del 70% del largo de la laguna se tiene que: 

Xf =(Lprom)(0.70(N 


+1) /((Bprom) /(N +1)) 


MampF MampF 


Luego, para el coeficiente de dispersión se 
tiene: 


df = y (21) 
-0.26118+0.25392(Xf) +1.0136(Xf) 


j) Coeficiente de reducción bacteriana. 


El coeficiente de reducción bacteriana (expre- 
sión (22) del método de diseño) es una constante 
cinética que depende de la temperatura, la cual 
se basa en la fórmula de Arrhenius (Marais, 
1974; Mara, Pearson, Oragui, Arridge, « Silva, 
2001; Banda, 2007): 
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K, =0.841(1.075) (22) 
k) Constante “a”: 
a,=[1+4(K, Ofdf) (23) 


1) Coliformes fecales en el efluente de la 
laguna facultativa: 


NE _ afec 920 


No  (l+afy (24) 


m) Coliformes fecales corregidos 


evaporación. Aplicar la fórmula (12). 


por 


n) Concentración de la DBO en el efluente de 
la laguna: 


ksa 


- (1.085) el 


Concentración de la demanda bioquímica de 
oxígeno en el efluente de la laguna: 


DBO'f 


DBO€f = 26 
E rOf+1 pa 
o) Eficiencia de remoción de la DBO: 
%= (DBO/¡f - DBOEf) x100 (27) 


DBO'f 


p) DBO corregida por evaporación. Aplicar la 
ecuación (11). 


Laguna de maduración (método flujo 
disperso) 


1. Tiempo de retención hidráulico. (O) es un 
método de prueba y error, donde se propone 
un tiempo de retención hidráulico y se 
calcula la concentración de los coliformes 
fecales en el efluente; esto, siguiendo el 
criterio de diseño de la laguna facultativa. 


2. Volumen de la laguna: 


Vm =(QD(Om) (28) 


3. Área de la laguna: 
(29) 


Para el resto del diseño se aplicaron las 
ecuaciones a partir del inciso f) hasta p), que se 
indica en la metodología de la laguna facultativa. 


Construcción del modelo matemático 


Los sistemas de tratamiento mediante lagunas 
de oxidación son sencillos tanto de construir 
como de operar, pero tienen la desventaja de 
requerir una gran cantidad de terreno para su 
construcción, en comparación con otros siste- 
mas más tecnificados, por lo que se propone 
un modelo matemático de minimización del 
área total ocupado por el sistema lagunar. Las 
restricciones y función objetivo deben estar re- 
lacionadas con las mismas variables cambiantes. 
Con el propósito de probar lo anterior, es nece- 
sario hacer las sustituciones correspondientes. 
El modelo propuesto, como ya se indicó, incluye 
una función objetivo, la cual se integra por la 
suma de áreas de las tres lagunas. Para las 
restricciones de la calidad del agua residual se 
consideró lo que indica la norma oficial mexi- 
cana NOM-001-SEMARNAT-1996. Los límites 
máximos permitidos de concentración son los 
siguientes: para coliformes fecales deben ser 
menores o iguales que 1 000 NMP/100 ml, y 
para la DBO, menores o iguales que 75 mg/]. 
En la figura 1 se indica el valor mínimo 

para el establecimiento de flujo disperso: 0.05. 
Este valor fue considerado tanto para la laguna 
facultativa como para la de maduración. Según 
Arceivala (1973), el tiempo de retención de un 
estanque facultativo varía de 10 a 110 días; 
mientras que para una laguna de maduración, 
el rango es mayor o igual que cinco días (WHO- 
EMRO, 1987). De esta manera se determinaron 
las restricciones. 

A continuación se muestra el modelo mate- 
mático en su forma reducida: 


Minimizar: 
A, =Aa+Af+ Am (30) 
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Sujeto a: 
DBO em <75 (31) 
Ne <1000 (32) 
d,=0.05 (33) 
d,, =0.05 (34) 
7<Of<110 (35) 
5<0m=<10 (36) 


Para llevar a cabo la demostración de las 
funciones del modelo matemático se iniciará con 
la restricción del número de coliformes fecales, 
por lo que es necesario sustituir en la ecuación 
(24) el número de coliformes en la entrada de la 
laguna facultativa por la ecuación (10), para así 
obtener la función (37): 


yin /2óo( Nc 
(1+4,) 1+k,,0, 

La ecuación (37) expresa el número de 
coliformes fecales en el efluente de la laguna 
facultativa; esta expresión se sustituye por el nú- 
mero de coliformes fecales también en la salida 
de la laguna de maduración en la fórmula (24) 
y resulta la expresión (38): 


de _%a, exp(1-a, /2d,) 4n¿exp(1-a,/2d/) 
en (l+a y (L+a,) 


m 


¿_N; (38) 
1+ K, O, 


La fórmula (38) se debe corregir por evapora- 
ción, multiplicándola por el cociente del caudal 
en la entrada de la laguna anaerobia entre el 
caudal en efluente de la laguna de maduración; 
de esta forma resulta la expresión (39): 


= 40 ,, exp(1-a2,, /2d )  An,exp(1-a,/2d,) 


N m 
ds (1+a, ) (1+a,) 
1 + Ko, O, 5 
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La fórmula (39) se puede reducir en exten- 
sión y resulta la función (40): 


d/(1-a,,)+d,, (1-a,) 


N 2 mty 1000 
= < 
em 2 2 ul 
Q., (1+a,,) *(1+a,) *(1+k,,0,) (40) 


lam 


16N,Q,,A,,4 0o| 


Los factores a, y a, equivalen a la fórmula 
(23) del método tradicional de diseño, aplicados 
a las lagunas facultativa y de maduración. Los 
factores d,y d,, equivalen a la expresión (21) del 
método de diseño aplicado a las lagunas faculta- 
tivas y de maduración. Al realizar la sustitución 
para llegar a la fórmula (40), se puede deducir 
que ésta depende del tiempo de retención de las 
lagunas facultativa y de maduración, así como 
del número de mamparas. La dispersión (d) 
depende de la relación largo ancho (X), que a 
su vez depende del número de mamparas y las 
dimensiones del estanque (Bprom y Lprom). Las 
dimensiones del estanque dependen del tiempo 
de retención; así se demuestra el vínculo entre la 
restricción y las variables de decisión. 

Es lo mismo para la restricción de los coli- 
formes fecales. Se deben llevar a cabo las susti- 
tuciones para la restricción de la DBO €. Para lo 
anterior se toma la expresión (26) y se sustituye 
la DBO en el influente (DBO/) por la fórmula de 
la DBO en el efluente de la laguna anaerobia y 
resulta la fórmula (41): 


100 
KfO; +1 


o, a] 


(41) 


DBO,, = 


La expresión (41) se sustituye en la fórmula 
(26) por la DBOI, para obtener la DBOe en el 
efluente de la laguna de maduración, se corrige 
por evaporación y resulta la fórmula (42), la cual 
representa la restricción: 


Do 
o Rola 75 (42) 
em (KfO, +11(KfO,, +1) 2 


La DBO? (expresión (42)) depende del tiem- 
po de retención en la laguna facultativa y de 
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maduración, variables de decisión del modelo 
matemático. 

Igual que las restricciones, la función objeti- 
vo (expresión (30)) depende de las variables de 
decisión. Lo anterior se puede demostrar si se 
sustituye en el área de la laguna facultativa y de 
maduración por el volumen dividido entre la 
profundidad. De esta forma resulta la expresión 
(43): 

Ad (43) 
Lo Za 


Se puede despejar el área de la ecuación (13) 
para la laguna anaerobia, a fin de demostrar que 
la función objetivo depende de la relación largo- 
ancho de la laguna anaerobia. De esta manera 
resulta la expresión (44): 


AE E 


t a” prom Z 


(44) 
f Z m 


Tomando como base la fórmula (28), se sus- 
tituye el volumen por el cociente del gasto en 
el influente, entre el tiempo de retención para 
las lagunas facultativas y de maduración, para 
obtener la fórmula (45): 


AX Bon. + Ly ¿Lo (45) 


a” prom 
0%; mm 


La fórmula (45) demuestra la relación entre 
las variables de decisión: relación largo-ancho 
en la laguna anaerobia (X), tiempo de retención 
en las lagunas facultativa y de maduración (O 
yO 

Se llevó a cabo un diseño tradicional para 


f 


,), con la función objetivo (4). 

un sistema de tres lagunas para una comunidad 
rural del municipio de Gómez Palacio, Duran- 
go. Con un gasto de diseño en la entrada de la 
laguna de 230 m? / día, una DBO, de 340 mg/l y 
coliformes fecales de 14 000 000 NMP/100 ml. 
Se consideró la temperatura promedio del mes 
más frío: 11.8 *C. Las profundidades conside- 
radas en las lagunas fueron 4.0 metros para la 
anaerobia, 1.5 para la facultativa y 1.0 metro 
para la de maduración o pulimiento. 


Resultados y discusión 


Los resultados del cálculo con el método tra- 
dicional de diseño se muestran en los cuadros 
1,2y3. 

De acuerdo con los resultados del cuadro 2, 
la dispersión (4) del estanque facultativo resultó 
de 0.1145. Este valor se encuentra en el rango 
de 0.05 y 8.0 (figura 1). Por lo tanto, el régimen 
hidráulico de la laguna es de flujo disperso (Pol- 
prasert £: Bhattarai, 1985; Aldana, Lloyd, Gu- 
ganesharajah, €: Bracho, 2005). Los coliformes 


Cuadro 1. Resultados del dimensionamiento de la laguna anaerobia con el método tradicional de diseño. 


Qia DBOia T Za ev Ni 
(Xa) Talud 
(m?/d) (mg/l) (eC) (2al4m) (mm/d) (NMP/100 ml) 
230 340 11.8 4 5 14 000 000 2 1 
Carga ; ; 
Carga a % DBO, A Tiempo medio 
a volumétrica id Volumen Área A eretenció DBOea Qea 
reánica removi e retención 
Bs da de diseño aa (m?/d) (m?) (Om) an (mg/1) (m?/d) 
a 
5 (gDBO,/m*-d) 
78.20 136.00 43.60 575.00 143.75 2.50 191.76 229.28 
N ta B D B D A 
ea : prom prom sup sup a 
Kt d 
kl (NMP/100 mb) | 800 (m) (m) (m) (m) (m) 
(mg/1) 
0.62 5 466 419 192.36 8.48 16.96 12.48 20.96 261.48 
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Cuadro 2. Resultados del dimensionamiento de la laguna facultativa con el método tradicional de diseño. 


Datos de entrada 
Qif Nea DBOjf 
(m?/d) (NMP/100 ml) (mg/l) mampf aña 
229.28 5 466 419 192.36 1 22.53 
Resultados 
Bprom Lprom Área 
Laguna Xf df Kb af P P promedio 
(m) (m) E 
(m?) 
Facultativa 8.4 0.1145 0.4648 2.40715 33.88 101.65 3 444.05 
PA Qef Nef DBOef T Bsup Lsup Af 
8 (m?/d) (NMP/100 ml) (mg/l) (20) (m) (m) (m?) 
Facultativa 209.98 10 592 41 11.8 36.88 104.65 3 859.64 


Cuadro 3. Resultados del dimensionamiento de la laguna de maduración con el método tradicional de diseño. 


Datos de entrada 
Qim Nef DBOim 
(m/d) (NMP/100 ml) (mg/D) amp alias 
209.98 10592 41 S 6.10 
Resultados 
Bprom Lprom Área 
Laguna Xm dm Kb am P P promedio 
(m) (m) 5 
(m?) 
Maduración 12.4963 0.0775 0.4648 1.37083 33.88 37.80 1 280.90 
aa Qem Nem DBOem T Bsup Lsup Am 
8 (m?/d) (NMP/100 ml) (mg/l) (eC) (m) (m) (m?) 
Maduración 202.46 980 20 11.8 35.88 39.80 1 428.27 


fecales resultaron todavía por encima de lo que 
indica la norma, por lo que es necesario agregar 
la laguna de maduración. La materia orgánica 
resultó por debajo del límite permitido. En el 
cuadro 3 se muestran los resultados de la laguna 
de maduración. 

De igual forma que la laguna facultativa, el 
estanque de maduración resultó con una dis- 
persión (4) para flujo disperso. Los coliformes 
fecales resultaron menor que 1 000 NMP/100 
ml, por lo que se cumple con la normatividad. 
La materia orgánica que se obtuvo también 
cumple con los valores permisibles de concen- 
tración (DOF, 1996). En la figura 2 se indican 
las dimensiones del diseño con el método 
tradicional. 
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Aplicación del modelo matemático 


Se aplicó el análisis matemático para la optimi- 
zación del área del sistema lagunar. Para llevar 
a cabo lo anterior, se utilizó el programa Solver 
de Excel. Como ya se indicó, éste considera el 
gradiente reducido generalizado (GRG). En el 
cuadro 4 se muestran los resultados de la laguna 
anaerobia optimizada. 

A partir de la metodología para el diseño 
de la laguna anaerobia, el porcentaje de la 
DBO, removida y la constante de decaimiento 
para los coliformes fecales (cuadro 4) dependen 
solamente de la temperatura (Conagua-IMTA, 
2007a, 2007b). Por lo tanto, la única manera de 
reducir el área del estanque es modificando la 
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20.96 m 


Efluente 


có 


12.48 m 


Maduración Facultativa 


35.88 m 36.88 m 


39.80m 104.65 m 


Figura 2. Dimensiones del sistema lagunar con el método tradicional de diseño. 


Cuadro 4. Resultados del diseño de la laguna anaerobia optimizada. 


. ; Ni 
sed EN e e a AS 1d) a (Xa) Talud 
230 340 11.8 4 5 14 000 000 1.005 1 
Carga a de %o BO, 5 Volumen Área en DBOea Qea 
E diseño senado (m/d) (m) retención (mg/D) (m/d) 
e (gDBO,/m'-d) (Oa) 
78.20 136.00 43.60 575.00 143.75 2.50 191.76 229.28 
Kta Nea a Bprom Lprom Bsup Lsup Aa 
(NMP/100 ml) (mg/D (m) (m) (m) (m) (m) 
0.62 5 466 419 192.36 11.96 12.02 15.96 16.02 255.67 
relación largo-ancho. El resultado obtenido por incrementó de 1 a 2. Al respecto, Muttamara, 2 
el modelo fue de un día. Según Senzia, Mayo, y Puetpaiboon (1996, 1997); Sperling, Cherni- E 
Mbwette, Katima y Jargensen (2002); Agun- charo, Soares y Zerbini (2002); Shilton y Ha- S 
wamba, Ochonogar y Ukpong (2003); Mara rrison (2003); Shilton y Mara (2005); Abbas et pa 
(2004); Abbas et al. (2006); Hamzeh y Ponce al. (2006); Banda (2007); Winfrey et al. (2010), y 3 
(2007), y Naddafi et al. (2009), la temperatura es Cortés-Martínez et al. (2014a; 2014b; 2015; 2016), E 
un factor importante en los sistemas de trata- publicaron que al utilizar mamparas en lagunas El 
miento de lagunas de estabilización y funcionan de estabilización se incrementa de manera im- De 
de mejor manera en climas tropicales. portante la eficiencia en la eliminación de conta- E 
En el cuadro 5 se muestran los resultados de minantes. Según Shilton y Mara (2005), y Abbas 5 
la laguna facultativa mejorada. et al. (2006), al utilizar mamparas en lagunas de E 
Se estableció la condición óptima con base estabilización se favorece un adecuado régimen > 
en las limitaciones planteadas en el modelo hidráulico, por lo que aumenta la remoción de 3 
matemático. contaminantes. Los resultados del presente i 
De acuerdo con los cuadros 2 y 5, el nú- estudio concuerdan con las afirmaciones de los 3 
mero de mamparas en la laguna facultativa se citados autores. S 
8 
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Cuadro 5. Resultados del diseño de la laguna facultativa optimizada. 


Datos de entrada 
Qif Nea DBO¡f Of (d) 
(m/d) (NMP/100 ml) (mg/l) mampf 
229.28 5 466 419 192.36 2 18.79 
Resultados 
Bprom Lprom Área 
Laguna Xf df Kb af pro pro promedio 
(m) (m) 
(m?) 
Facultativa 18.9 0.0516 0.4648 1.67359 30.94 92.82 2872.11 
auna Qef Nef DBOef T Bsup Lsup Af 
5 (m¿/d) (NMP/100 ml) (mg/l) (e0) (m) (m) (m?) 
Facultativa 213.02 8 016 47 11.8 33.94 95.82 3 252.41 


La concentración de los coliformes fecales 
obtenidos resultó por encima de los límites 
máximos de concentración que se indican en 
la norma oficial mexicana. Con esto se justifica 
agregar una laguna de maduración o pulimien- 
to. La DBO, resultó por debajo del límite de la 
norma. En el cuadro 6 se muestran los resul- 
tados de la laguna de maduración optimizada. 

Con base en el cuadro 6, el flujo resultó 
disperso 0.0501. Según Arceivala (1981) y 
Mara (2004), algunos factores que favorecen la 
implementación del flujo disperso son la pro- 
fundidad, la ubicación de entradas y salidas del 


flujo, y el número de mamparas o deflectores. 
Como se observa en los cuadros 3 y 6, se indica 
el número de mamparas, que se incrementaron 
de 3 a 4; pero en ambos análisis se mantiene 
el flujo disperso. En la figura 3 se muestran las 
dimensiones del sistema lagunar, calculadas con 
los criterios de optimización. 

En el cuadro 7 se muestran los resultados de 
las áreas superficiales, así como los tiempos de 
retención calculados tanto con la metodología 
tradicional como con el modelo matemático 
de optimización. En la laguna anaerobia, como 
ya se dijo, al modificar la relación largo-ancho 


Cuadro 6. Resultados del diseño de la laguna de maduración optimizada. 


Datos de entrada 
Qim Nef DBOim 
(m/d) (NMP/100 ml) (mg/l) mampm emo) 
213.02 8016 47 4 5.0 
Resultados 
Bprom Lprom Área 
Laguna Xm dm Kb am P P promedio 
(m) (m) s 
(m?) 
Maduración 19.4691 0.0501 0.4648 1.2105 30.94 34.42 1 065.10 
a Qem Nem DBOem T Bsup Lsup Am 
8 (m*/d) (NMP/100 ml) (mg/l) (0) (m) (m) (m?) 
Maduración 206.67 1 000 20 11.8 32.94 36.42 1 199.82 
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Anaerobia 


Efluente 


E —> 
Noel 
2 7] : ES 
E tuene Facultativa ES 
16.02 m- 95.82 m 36.42 m 
Figura 3. Dimensiones del sistema lagunar con el modelo de optimización. 
Cuadro 7. Resumen de áreas y tiempo de retención en el diseño y optimización del sistema lagunar. 
A 2% e Diferencia tiempo E o 
Lagunas Método tradicional Modelo matemático e Diferencia área 
de retención 
Tiempo de Tiempo de 
Área (m3) retención | Área (m?) retención (d) % (m2 % 
(d) (d) 
Anaerobia 261.48 2 255.67 25 0 0 5.81 LD 
Facultativa 3 859.64 22.53 3 252.41 18.79 3.74 16.60 607.23 15.73 
Maduración 1 428.27 6.10 1 199.82 5 11 18.03 228.45 15.99 
ds 5 549.39 31.13 4707.9 26.29 4.84 841.49 15.16 
porcentajes 


se obtuvo una reducción de 2.22% de área; 
mientras que en la laguna facultativa resultó 
menor tiempo de retención y área con el modelo 
matemático, 16.60 y 15.73%, respectivamente. 
Por último, en la laguna de maduración, el 
área se redujo 15.99%. En total, la reducción del 
área y tiempo de retención con el modelo de 
optimización resultó de 4.84 días y 841.49 me- 
tros cuadrados, que representa 15.16%. Según 
Rolim (2000) y Conagua-IMTA (2007a, 2007b), la 
principal desventaja de los sistemas lagunares 
es la gran área de terreno que requieren para su 
construcción. 

La novedad, en comparación con las pu- 
blicaciones de los autores Olukanni y Ducoste 
(2011), y Cortés-Martínez et al. (2013, 2014a, 
2014b, 2015, 2016), es un modelo matemático 
integrado por tres lagunas, con restricciones de 
flujo disperso; el objetivo de la optimización es 
el área que ocupa el sistema. 


Análisis de sensibilidad 


De acuerdo con Saltelli et al. (2008), el estudio 
consiste en realizar una modificación de los va- 
lores de entrada del modelo (X, O, Om, N 
) en incrementos de + 2%, con el propósito 


mampf” 


mampm 


de observar los cambios en los resultados con 
respecto a los valores optimizados (0% de va- 
riación). La figura 4 muestra los resultados del 
análisis de sensibilidad. 

En la figura 4, la coordenada (0.00, 1.00) 
indica los valores determinados por el modelo 
matemático. A partir de este punto se observa 
que cuando los parámetros de diseño se incre- 
mentan de 0.00 a 0.06, la demanda bioquímica 
de oxígeno y coliformes fecales decrecen, a 
diferencia de A, que aumenta. Lo anterior se 
puede interpretar de la siguiente manera: a 
mayor tiempo de retención mayor área. Cuando 
los parámetros decrecen, la tendencia es inversa: 
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Incremento con respecto a la optimización 


-0.06 -0.04 -0.02 -0.00 


—— Ne 
-—a- DBOe 
4 AT 


mb 1.05 


0.02 0.04 0.06 


Variación de los parámetros de diseño 


Figura 4. Análisis de sensibilidad del modelo matemático. 


a menor tiempo de retención y número de mam- 
paras mayor concentración de contaminantes 
(DBOe y Ne). 


Conclusiones 


De acuerdo con los objetivos planteados con la 
aplicación del modelo matemático propuesto 
fue posible reducir el área en un sistema inte- 
grado por tres lagunas. El aporte del presente 
trabajo fue incluir la restricción para flujo 
disperso en el modelo matemático. Se obtu- 
vieron mejores resultados: menor tiempo de 
retención y área. La diferencia con respecto a 
los resultados con la metodología tradicional 
es de consideración. 

Es importante indicar que el criterio de 
optimización para la laguna anaerobia fue 
diferente: se consideró la relación largo-ancho. 
Para la laguna facultativa y de maduración se 
aplicaron las restricciones incluidas en el mode- 
lo matemático. 

El análisis de sensibilidad usado en el mo- 
delo matemático resultó favorable al cambio 
de las variables. El criterio y análisis para la 
aplicación del modelo matemático propuesto se 
puede emplear en diferentes condiciones, pero 
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se necesitan cambiar los siguientes parámetros: 
temperatura promedio del mes más frío y eva- 
poración. 
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